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ABSTRACT 

The  development  of  a  pho t ogr amme t ry  software  package 
requires  considerable  study.  Computer  design  and  capacity, 
available  phot ogr ammet rlc  instruments,  analytical  methods 
versus  semi-analytical  methods,  simplicity  versus  compre¬ 
hensiveness,  user  competence,  and  the  desired  final  products 
are  but  a  few  of  the  elements  to  be  weighed.  An  example 
system  is  presented  for  consideration  in  this  paper.  This 
system  is  being  implemented  throughout  Latin  America  by  the 
Inter  American  Geodetic  Survey.  The  package  consists  of 
both  an  analytical  and  a  semi-analytical  adjustment  program 
and  the  accompanying  programs  which  tie  the  systems  together 
and  allow  the  user  to  go  from  relative  orientation  or  com¬ 
parator  measurements  to  plotter  element  settings  via  computer. 
Included  in  the  discussion  is  control  selection  and  distri¬ 
bution,  pass  point  selection,  program  execution,  and  analysis 
of  results.  The  sample  software  package  was  assembled 
because  it  is  simplC^MiMr  and  yet  versatile  enough  to  be 
adapted  to  a  variety  of  computers  and  phot ogr ammet ric  equip¬ 
ment.  It  is  presented  as  a  successful  example  of  a  system 
design,  and  it  is  hoped  that  this  discussion  and  example 
will  generate  new  ideas  in  system  development. 

INTRODUCTION 


Photogr ammet ry  may  be  defined  simply  as  a  method  of  extract¬ 
ing  information  about  an  object  by  studying  its  image. 
Photographs  are  the  most  familiar  types  of  images,  but 
devices  sensitive  to  other  parts  of  the  electromagnetic 
spectrum  may  also  be  used.  Regardless  of  the  type  of  sen¬ 
sor,  the  image  must  be  analyzed  before  the  correct  informa¬ 
tion  is  obtained.  The  desired  results  could  be  interpretive 
and  require  the  art  of  pho t o -in t e rp r e t at  ion  to  identify  and 
evaluate  the  image.  However,  the  desired  results  could  be 
quantitative  and  require  measurements  and  mathematics  to 
arrive  at  the  desired  goal. 
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The  advent;  of  the  computer  and  its  adoption  to  photogram- 
metric  problems  has  greatly  enhanced  this  science.  Even 
in  photo-interpretation  the  computer  is  used  to  recognize 
patterns  in  even  the  nonvisibie  region  of  the  spectrum. 
Quantitative  methods  have  benefited  immensely  from  com¬ 
puter  utilization.  It  is  now  possible  to  solve  thousands 
of  simultaneous  equations  in  only  a  short  time. 

The  computer  software  discussed  in  this  paper  was  designed 
to  help  obtain  quantitative  results  from  aerial  photography. 
The  results  are  used  to  obtain  dimensions  and  precise 
positions.  The  desire  for  increased  accuracy  has  led  to 
the  development  of  precise  equipment,  systematic  procedures, 
and  rigorous  computer  adjustments.  The  most  common  appli¬ 
cation  of  the  resultant  information  is  in  the  preparation 
of  topographic  maps  and  charts  of  the  earth's  surface. 

Almost  all  of  the  national  mapping  programs  in  the  Western 
Hemisphere  are  performed  by  pho t ogr amme t r y . 

The  goal  of  the  software  system  developed  by  IAGS  is  to 
encompass  the  needs  of  all  the  associate  agencies.  It  is 
not  burdened  with  copyright  laws  and  therefore  allows 
everyone  free  access.  It  maintains  simplicity  and  yet 
allows'us  to  obtain  meaningful  results  from  a  wide  variety 
of  equipment.  Technical  advancement  in  computer  hardware, 
photogrammetric  instruments,  and  mathematical  formulation 
will  not  allow  any  software  system  to  rest.  Our  programs 
are  under  constant  change,  and  we  are  always  looking  for 
something  better.  The  package,  however,  has  been  installed 
and  is  being  successfully  used  in  the  United  States  and  in 
several  Latin  American  countries. 

SYSTEM  DESIGN 

The  design  of  any  photogrammetric  software  package  must 
begin  with  the  adjustment  program.  It  is  the  largest  and 
most  complex  part  of  the  system,  and  the  other  programs 
are  built  around  it.  The  wide  variety  of  instrumentation, 
technology,  and  requirements  found  in  the  hemisphere  pre¬ 
clude  the  use  of  a  single  program  to  accomplish  this  task. 
Thus  both  semi-analytical  and  fully  analytical  approaches 
were  taken.  Both  applications  produce  accuracy  and  effi¬ 
ciency  sufficient  to  fulfill  our  mapping  requirements. 

The  purpose  of  the  adjustment  program  is  to  enhance  the 
ground  control  network  that  currently  exists  through  field 
surveying,  doppler  positioning,  or  other  techniques.  Each 
model  used  in  a  stereoplotter  requires  at  least  two  hori¬ 
zontal  control  points  for  scaling  and  three  vertical  control 
points  for  leveling.  It  would  not  be  economical  to  estab¬ 
lish  all  this  control  by  ground  methods.  Adjustment  programs 
must  also  provide  a  method  of  detecting  blunders  in  the 
observations,  evenly  distribute  small  systematic  errors, 
and  determine  estimates  of  precision  for  the  results. 

After  the  adjustment  programs  were  selected,  a  series  of 
auxiliary  programs  were  acquired  and  written  to  blend  the 
two  main  programs  into  a  versatile  smooth-flowing  system. 


These  programs  allow  us  to  check  observations  for  blunders, 
cheek  data  dimensions  so  as  not  to  exceed  the  adjustment 
parameters,  produce  independent  model  coordinators  from 
two  dimensional  comparator  measurements,  mathematically 
compute  projection  centers,  determine  stereo plotter  orien¬ 
tation  settings  after  the  adjustment,  and  provide  other 
useful  information  as  will  be  discussed  later  in  this 
paper. 

A  depiction  of  the  overall  system  flow  is  shown  in  Figure  1. 
Raw  data  may  enter  the  system  from  three  different  sources, 
be  channeled  in  various  directions  and  adjusted  as  desired. 

THE  SIMBA  ADJUSTMENT  SYSTEM 

The  Simultaneous  Block  Adjustment  of  Models  (SIMBA)  program 
was  written  by  Randle  W .  Olsen  of  the  U.S.  Ceological  Survey 
in  1973.  It  has  been  adapted  for  use  by  the  U.S.  Forest 
Service  and  DMA  IAGS.  It  is  written  in  FORTRAN  IV  and  is 
used  to  simultaneously  adjust  independent  model  units  to 
each  other  and  to  ground  control.  The  three-dimensional 
data  may  be  generated  on  a  stereoplotter  or  on  a  comparator 
providing  the  CRAST  (Coordinate  Re f inement -Ana ly t i cal  Strip 
Triangulation)  or  a  similar  program  is  used  to  derive  the 
projection  centers  and  the  third  dimension.  The  models  are 
adjusted  by  an  interactive  series  of  planimetry-height  solu¬ 
tions.  Four-parameter  linear  transformations  are  used  in 
the  horizontal  solution  and  three-parameter  linear  trans- 
fromations  are  used  in  the  vertical  solution.  This  separation 
economizes  computer  time  and  storage  without  a  significant 
compromise  in  accuracy  as  compared  to  more  traditional  seven- 
paramater  adjustments. 

Program  features  include: 

1.  Input  arranged  in  model  units  with  a  separate  control 
deck  . 

2.  Internal  sorting  of  control  and  tie  points. 

3.  In-core  solution  of  banded  normal  equations  to  minimize 
computer  time. 

4.  Separate  weighting  for  ground  control  and  model  ties. 

5.  Printed  output  arranged  in  model  units  with  corres¬ 
ponding  residuals. 

6.  A  listing  of  test  point  coordinates  and  residuals  which 
were  withheld  from  the  solution. 

7.  A  root  mean  square  (RMS)  error  summary  of  held  control, 
tie  points,  and  unheld  control. 

8.  Absolute  orientation  values  to  include  common  tilt, 
common  tip,  the  airbase,  and  exposure  station  elevation 
difference  for  each  model. 


Special,  attention  should  he  given  to  the  weighting  system 
as  it  affords  us  the  opportunity  to  separate  errors  result¬ 
ing  from  poor  control  from  errors  due  to  relative  orientation 
or  had  model  tie  points, 
of  the  weight  ratio  is: 

W  =  - n - j- 

Sm -  +  S  c  ^ 

Where  Sm  =  estimated  unit  coordinate  standard  deviation 
of  the  model  points 

Sc  =  estimated  standard  deviation  of  each  control 
point 

One  can  easily  see  that  the  value  of  ”W"  will  range  from 
zero  to  two.  A  value  of  one  will  allow  equal  weight  to 
be  placed  on  the  control  value  and  the  model  ties.  A 
value  close  to  two  will  constrain  the  control  and  allow 
the  models  freedom  to  move.  Likewise  a  value  close  to 
zero  will  constrain  the  models  and  allow  the  control  to 
shift  . 

A  minimum  of  two  horizontal  control  points  are  required 
for  t  h-e  planimetric  adjustment,  and  three  vertical  control 
points  are  required  for  the  height  adjustment.  If  no 
additional  control  is  used,  there  will  be  an  absolute 
solution  and  no  errors  will  be  propagated  by  a  least 
squares  adjustment.  Thus  if  only  minimum  control  is 
used  and  constrained  while  the  model  ties  are  given  free¬ 
dom  to  move,  the  propagated  error  will  result  from  poor 
model  fit,  and  the  tie  point  residuals  will  be  exaggerated. 
Actual  errors  present  in  the  constrained  control  are  un¬ 
important  because  they  are  only  being  used  to  scale  and 
level  the  block  solution,  and  we  are  only  checking  tie 
point  residuals. 

Once  a  good  model  tie  solution  is  achieved,  the  full  control 
network  is  used,  but  the  weight  is  changed  to  constrain  the 
models  and  give  freedom  of  movement  to  the  control.  The 
control  residuals  thus  become  exaggerated,  and  control 
values  incompatible  with  the  relative  model  solution  can 
quickly  be  noted.  After  all  blunders  have  been  eliminated, 
equal  weighting  is  generally  used  in  the  final  solution 
unless  extenuating  circumstances  dictate  otherwise. 

The  computer  memory  necessary  to  execute  this  adjustment  is 
largely  dependent  upon  the  size  of  the  normal  equation 
coefficient  matrix.  If  the  independent  models  are  properly 
aligned  in  the  program,  this  matrix  becomes  banded  along 
the  diagonal  leaving  only  zeros  outside  the  band.  To  con¬ 
serve  storage,  only  the  elements  within  the  band  above  the 
main  diagonal  are  stored.  They  are  arranged  in  a  rectangular 
array  called  "Q"  in  the  program  and  are  dimensioned  4  x  (the 
bandwidth)  by  4  x  (the  number  of  models).  Correct  sequencing 
of  the  models  is  therefore  critical  to  stay  within  the 
allotted  bandwidth.  In  a  rectangular  block,  the  band¬ 
width  can  roughly  equal  the  number  of  models  in  the  longest 
strip  plus  two.  The  normal  equations  are  solved  by  the 
Gaussian  elimination  procedure  using  routines  suited  for 
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symmetric  arul  baruied  equations.  The  banded  routine  con¬ 
serves  computer  core  and  significantly  reduces  execution 
time. 

It  is  always  difficult  to  recommend  a  control  point  pattern 
because  there  arc  so  many  variables,  the  most  important  of 
which  is  quality.  Also,  it  has  been  ray  belief  that  a  good 
pho  t  og  r  amme  t  r  is  t  will  attempt  to  secure  as  much  control 
information  as  possible,  and  his  only  limiting  factor  will 
bo  tlie  man  in  charge  of  expenditures.  At  the  risk  of 
violating  this  rule,  a  control  pattern  which  should  satisfy 
most  needs  is  recommended. 

It  is  essential  to  have  horizontal  control  in  the  extreme 
corners  of  the  block.  Any  deviation  from  this  will  cause 
warping  in  the  corners  becoming  progressively  worse  as  the 
control  point  strays  from  its  corner  position.  Other  than 
the  corners,  SIMBA  seems  to  hold  horizontal  position  rather 
well.  Additional  horizontal  control  should  be  spaced 
around  the  edges  of  the  block  at  six  to  eight  model  inter¬ 
vals. 

Vertical  control  is  a  different  matter.  The  adjustment 
acts  4,s  a  blanket  tacked  down  at  the  corners  with  a  fan 
blowing  air  under  it.  It  bulges  in  the  middle,  and  when 
we  tack  down  the  middle,  we  develop  two  smaller  bulges  on 
each  side.  This  progression  continues  geometrically  until 
the  bulges  are  reduced  in  size  to  an  acceptable  level.  The 
general  recommendation  is  to  have  a  vertical  control  pattern 
at  three  to  four  photo  intervals  throughout  the  block.  A 
sample  control  point  pattern  is  shown  in  Figure  2. 

Pass  and  model  point  selection  is  a  matter  of  personal 
custom,  but  the  author  has  developed  a  system  which  lends 
itself  well  to  both  the  adjustment  mathematics  and  production 
flow.  Starting  with  one  of  the  border  strips  in  the  block, 
pick  four  points  in  a  line  in  all  the  trilap  areas  and  at 
the  two  ends.  Insure  that  the  two  points  closest  to  the 
adjacent  strip  fall  in  the  sidelap  area  and  are  trans¬ 
ferred.  Transfer  these  tie  points  to  the  second  strip 
and  send  the  first  strip  to  the  comparator  or  stereoplotter 
to  be  measured.  Pick  points  on  the  second  strip  as  before. 

If  transferred  points  from  the  first  strip  happen  to  fall 
in  the  trilap  or  necessary  end  positions  of  the  second 
strip,  use  them.  If  not,  continue  point  selection  as  on 
strip  one.  It  is  not  necessary  to  transfer  any  of  the 
newly  selected  points  back  to  the  first  strip.  Continue 
in  this  manner  until  all  points  are  selected.  This  method 
will  produce  from  eight  to  ten  pass  points  per  model  and 
is  illustrated  in  Figure  3. 

Three  versions  of  SIMBA  are  currently  available  at  DMA  IAGS 
Headquarters  in  San  Antonio,  Texas.  These  versions  are  the 
standards  which  are  then  adapted  to  the  computer  and  desires 
of  individual  users.  The  first  version  is  the  most  common 
and  has  parameters  of  240  models,  2200  unique  points,  2500 
total  points,  150  control  points,  and  a  bandwidth  of  25. 

It  uses  526K  bytes  of  storage  on  an  IBM  370  computer.  The 
second  version  of  SIMBA  breaks  the  program  down  into  three 
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purrs  wlii  i'll  -ire  ox  ecu  tod  separately.  The  first  part  reads 
■ho  in:>ut  data,  organizes  it,  and  builds  a  file.  The 
-.ei-oiid  part  reads  that  file,  performs  the  necessary  mathe¬ 
matics,  and  writes  the  massaged  data  back  onto  the  file. 

The  third  part  reads  the  file  and  prints  out  the  desired 
results.  Tii  is  three  part  version  is  a  computer  core  saver 
but  sacrifices  time  and  simplicity.  Dimensioned  similar 
to  tiie  original  program,  the  three  part  version  will  save 
about  75K  bytes  of  storage.  The  third  version  we  have  is 
another  three  part  program  with  the  first  and  tiiird  parts 
remaining  as  before.  The  difference  in  part  two  is  that 
a  direct  access  disk  file  is  used  to  write  and  then  read 
each  row  and  column  of  the  normal  equation  coefficient 
matrix.  All  of  this  reading  and  writing  takes  time  but 
an  additional  reduction  of  125K  bytes  of  computer  memory 
is  realized.  This  version  is  impractical  unless  the  computer 
memory  is  limited  and  computer  time  is  of  no  concern. 

A  sample  SIMBA  output  is  shown  in  Appendix  A. 

The  SIMBA  adjustment  requires  the  use  of  perspective  centers 
to  strengthen  the  vertical  solution.  If  a  stereoplotter  is 
used  for  relative  orientation,  the  coordinates  for  the  per¬ 
spective  centers  may  be  computed  using  either  one  of  two 
available  programs.  The  ideal  solution  for  determining 
the  camera  perspective  center  is  by  space  resection  using 
a  reseau  grid.  The  IAGS  version  of  John  McLaurin’s  Per¬ 
spective  Center  Determination  Program  (U.  S.  Geological 
Survey,  1969)  does  precisely  this.  A  grid  of  known  pre¬ 
cision,  usually  engraved  on  a  glass  plate,  is  projected 
through  a  stereoplotter  projector.  The  coordinates  of 
f  r-i  three  to  fifty  grid  Intersections  are  then  measured 
in  model  space.  The  bundle  of  rays  passing  through  the 
grid  intersections  on  the  calibrated  plate  and  the  bundle 
extending  to  the  same  measured  intersections  in  model  space 
originate  from  the  same  point  —  the  perspective  center.  If 
the  latter  bundle  of  rays  is  fitted  to  the  other  bundle  in 
a  least  squares  space-resection  adjustment,  accurate  coordi¬ 
nates  for  the  perspective  center  can  be  determined.  Printed 
output  includes: 

1.  Title. 

2.  Input  data  to  include  point  numbers,  number  of  read¬ 
ings  on  each  point  and  grid  coordinates. 

3.  Number  of  points  used  and  computed  principle  distance. 

4.  Mean  projected  grid  coordinates  and  standard  devia- 
t ions  . 

5.  Parameter  values  after  each  iteration. 

6.  Residuals  of  the  projected  points  in  their  coordinate 
sys  t  era . 

7.  Variance-covariance  matrix. 


8.  Standard  errors. 


Although  this  program  is  quite  accurate  and  produces  a 
wealth  oi  inioraat  ion,  it  in  time  con  sum  in/,  to  use.  This, 
coupled  with  the  unavailability  of  calibrated  grid  plates 
and  stereoplotter  variety  encouraged  IAC.S  to  develop  a 
simplified  version.  For  instruments  s u  c h  as  t  h  e  W  i  1 d  A -  7  , 
A-8,  A-9  and  A-10,  t  li  e  perspective  centers  do  not  vary  with 
the  relative  orientation  of  the  model  since  the  axis  of 
rotation,  for  ones; a  and  phi  intersect  the  projection  cardan. 
Tiie  perspective  centers  change  only  if  the  instrument  base 
component  is  changed.  Thus  tile  coordinate  determinations 
need  be  made  only  once  per  strip  or,  in  special  cases,  per 
block.  However,  for  optical  instruments  such  as  the  Zeiss 
C-.H  and  Wild  B-S,  the  projection  center  varies  with  the 
relative  orientation  of  each  model  since  the  axes  of  rota¬ 
tion  for  omega  and  phi  do  not  intersect  the  projection 
cardan.  If  the  perspective  center  is  to  be  determined 
after  each  model  relative  orientation,  a  quick,  method  is 
needed . 

Instead  of  a  reseau  grid,  the  engraved  lines  and  crosses 
on  the  plate  carriers  are  measured  and  substituted  for  a 
grid.  These  intersections  can  be  quickly  measured  after 
relative  orientation  of  each  independent  model.  The 
instrument  z  and. the  numbering  sequence  must  remain  con¬ 
stant  for  each  set  of  recordings.  The  program  is  designed 
to  compute  perspective  centers  of  the  left  and  right 
cameras  alternately.  A  maximum  of  ten  points  may  be  used, 
but  six  points  are  customary.  The  program  output  is 
limited  to  the  y,  and  z  coordinates  of  each  projection 

center,  out  this  is  sufficient  input  for  the  S IMB A  adjust¬ 
ments.  An  example  program  execution  is  shown  in  Appendix  B. 

One's  first  experience  with  a  large  adjustment  program  can 
be  frustrating  if  two  unique  points  have  the  same  point 
number.  The  block  tends  to  fold  in  upon  itself  causing 
large  residuals  and  an  unsatisfactory  solution.  Another 
frequent  error  is  the  exceeding  of  the  bandwidth  parameter. 
This  is  a  result  of  long  flight  lines  or  incorrect  model 
sequencing  in  the  program.  SIMBA  is  unpredictable  when 
this  occurs  and  may  diverge  instead  of  reaching  a  solution. 
These  and  various  other  errors  present  in  the  data  can  be 
difficult  to  locate.  Also,  it  is  expensive  to  run  large 
adjustments  if  errors  are  inherent  in  the  data.  For  these 
reasons  the  author  wrote  the  SIMBA  SORT  program  (SISRT). 

It  uses  the  exact  same  data  deck  (or  file)  as  SIMBA  and  is 
quick  and  inexpensive  to  execute.  The  main  program  function 
is  to  do  a  point  sort  and  list  the  location  by  strip  and 
model  of  each  unique  point.  It  also  computes  and  lists  the 
bandwidth  of  each  point.  This  allows  the  user  to  insure 
the  SIMBA  bandwidth  is  not  violated,  or  if  it  is,  which 
point  is  the  culprit.  A  complete  data  deck/file  listing 
is  provided.  Sample  program  output  may  be  found  in  Appen¬ 
dix  C . 

The  Numerical  Model  Orientation  Program  (ORELM)  was  written 
by  Randle  Olsen  of  the  U.  S.  Ceological  Survey  in  1976. 

The  purpose  of  the  program  is  to  use  the  omega  and  phi 
elements  from  relative  orientation,  the  desired  model 
scale,  and  the  common  omega,  phi,  airbase,  and  bz  for 
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the  U  a  l l  e  (  MX- S3)  ,  Kelsh,  W  i  1 .1  3-.J,  Korn  i’C-2  and  the 
•  •  »a  .»-■!,>  p  i’d  s  .  no  iii'  c>t  rront  i  v  wo  rk  in-,  to  incorporate 
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r  e  i;  i:  ;  r  ••  graduated  dial;  or  auxiliary  a  u  i  p  -i  •  ■  n  t  to  u  t  i  !  i  ’  • 
t  h  i  a  data.  The  use  o  l  the  program  e  I  i  ;n  i  a  a  r.  a  tie-  need  to 
reorient  ,  scale,  and  level.  1>  .*  t'oro  plot'  in;;.  Sample  output 

he  found  in  Append  i  ::  11. 

thi:  giant  au.ii*:;tment  systkm 

In  Cliia  analytical  ph  o  t  op,  r  arcane  L  r  i  c  system,  the  physical 
reconstruction  of  model  ;  e  o  m  e  t  r  y  is  replaced  with  mathe¬ 
matics.  T'ae  mathematical  model  is  constructed  to  represent 
the  relationship  e  t  w e  e n  the  photographic  images  and  the 
object  space.  The  images  are  measured  on  a  comparator  and 
numerical  methods  are  used  to  produce  the  position  and 
orientation  of  the  exposure  station  and  the  image  coordi¬ 
nates  in  object  space. 

Analytical  and  semi-analytical  methods  are  both  capable  of 
producing  accurate  results.  Therefore,  accuracy  is  not  the 
justification  for  the  extra  expense  of  a  comparator  and  the 
analytical  programs.  However,  there  are  inherent  advantages 
in  an  analytical  system.  First  of  all,  it  is  much  easier 
t  -  bridge  water  or  snevj  fields.  One  does  not  have  to  worry 
about  trying  tc  remove  paralax  from  a  model  that  is  half 
water.  Secondly,  a  large  variety  of  photography  and  cameras 
can  be  accommodated .  Panoramic  photography,  for  example, 
dees  not  have  a  simple  central  projection  and  special  cam¬ 
eras  on  high  flying  aircraft  may  have  a  very  long  focal 
length  which  is  impossible  to  duplicate  on  a  stereo  com¬ 
parator.  Thirdly,  analytical  methods  permit  the  input  of 
auxiliary  sensors.  Such  sensors  include  stabilization 
systems,  profile  recorders,  altimeters  and  astronomical 
observations.  This  additional  information  may  be  incorpor¬ 
ated  directly  into  the  solution  and  weighed  according  to 
reliability.  It  should  be  noted  that  some  semi-analytical 
programs  now  have  this  ability. 

Analytical  photogrammetry  offers  great  potential.  It  can 
incorporate  various  distortions  and  deformations  which  are 
impossible  to  duplicate  in  stereoplottcrs.  The  proponents 
of  analytical  methods  claim  they  should  be  able  to  obtain 
more  accurate  results,  but  in  practice  this  has  not  proven 
to  be  true.  The  complex  mathematics  and  equipment  costs 
have  frightened  away  many  potential  practitioners. I ACS  does 
not  recommend  one  method  over  the  other,  but  instead  has 
incorporated  both  methods  into  the  software  package. 

When  pho to gramme t r ic  images  are  measured  on  a  monoscopic 
or  stereoscopic  comparator,  two  dimensional  coordinates 
are  produced.  The  CRAST  program  refines  these  coordinates 
and  produces  either  GIANT  input  or  independent  model  coordi¬ 
nates  for  input  into  SI MBA.  Mathematical  refinement  of  the 
comparator  measurements  is  necessary  because  the  model  space 
was  not  physically  reconstructed.  It  is  mathematically 


reconstructed,  so  we  need  t  o  knew  additional  Information 
to  accurately  duplicate  what  took  place  when  the  photos 
were  taken.  The  first  three  items  on  the  following  list 
are  required  for  program  execution.  The  remaining  items 
will  be  u  s  e  ii  if  t  ti  e  y  are  known  and  input. 

1.  Camera  Focal  Length 

2.  Approximate  Flying  Height 

i  .  Sehut  Refraction  Coefficient 

♦  .  Comparator  Scale  Corrections  in  x  and  y 

5.  Comparator  Non- Orthogonality 

6.  Initial  Airbase  at  Photoscale 

7.  Radial  Lens  Distortion  Corrections 

8.  Principle  Point  Variation  from  Photo  Center 

9.  -Calibrated  Fiducial  Coordinates 

Depending  on  the  information  input  and  the  desires  of  the 
user,  the  CRAST  program  will  correct  the  input  coordinates 
for  the  following  distortions: 


1  .  Film  Distort  ion 

2.  Earth  Curvature 

3.  Light  Refraction 

4.  Comparator  Calibration  Corrections 

5.  Radial  Lens  Distortion 

If  the  GIANT  option  is  used,  CRAST  output  will  consist  of 
a  fit  of  fiducials  to  camera  data,  the  refined  plate  coordi¬ 
nate  listing  and  punched  card  output  in  the  GIANT  format  or 
an  optional  format.  If  the  SIMBA  option  is  used,  output 
will  consist  of  independent  model  coordinates,  y-parallax, 
model  tie  residuals  and  card  output  in  the  SIMBA  or  an 
optional  format. 

The  reason  for  the  use  of  this  particular  program  is  its 
adaptability  to  the  lack  of  information  which  frequently 
is  found  to  be  the  case  in  Latin  America.  For  example, 
the  program  will  accept  the  absence  of  calibrated  fiducials 
and  even  the  absence  of  fiducial  observations.  CRAST  was 
originally  written  by  Randle  Olsen  in  1973  and  has  been 
revised  several  times,  including  a  revision  by  IAGS  in 
1981.  Sample  program  output  can  be  found  in  Appendix  E. 


The  G  e  n  o  r  .1  l  Intesratoii  An  .1 1  y  t  i  <■  1 1  T  r  i  an;;  11  L  .1 1  i  on  Program 
(GIANT)  is  (losi  ;ni’d  to  perform  a  l  oast  squares  .id  j  ustnont 
ol  .1  block  of  frame  photu);r:iph:i  .  It  was  written  by  At  el  A. 

E  las  sal  in  1 9  7b  and  has  been  adapted  for  use  by  the  II.  S. 
Geological  Survey,  the  G .  S.  Forest  Service  and  I  AGS. 

Manv  variations  of  tiie  program  now  exist  and  at  this 
writing,,  it  is  currently  undergoing  yet  another  revision 
at  I AGS. 

GIANT  will  analytically  solve  for  the  ground  coordinates  of 
image  points  measured  on  two  or  more  photographs.  It  also 
solves  for  each  camera  station  position  and  orientation. 

Only  uncorrelated  observations  are  acceptable  and  may  be 
individually  weighed  to  reflect  known  precision.  This  system 
allows  the  mixture  ‘of  horizontal  and  vertical  control  of 
varying  accuracies.  It  also  allows  the  use  of  known  camera 
station  parameters.  The  program  propagates  error  estimates 
through  the  solution,  computes  the  a  posteriori  estimate  of 
the  variance  of  unit  weight  and  has  an  option  for  listing 
the  variance-covariance  matrix  and  standard  deviation  of 
the  output  parameters.  These  indicators  are  useful  esti¬ 
mates  of  the  solution  accuracy. 

An  initial  approximation  is  required  of  each  camera  station 
position  and  orientation  and  of  each  unique  image  point. 

The  approximations  of  the  image  points  are  computed  inter¬ 
nally  and  gross  approximations  are  acceptable  for  the  camera 
parameters . 

Program  dimensions  require  that  no  strip  exceeds  20  photo¬ 
graphs  or  that  there  are  fewer  than  ten  strips  in  the 
adjustment.  IAGS  currently  has  two  versions  of  GIANT  on 
file  with  the  difference  being  in  the  dimensions.  The 
principle  version  has  the  following  parameters: 

1.  460  photographs 

2.  450  ground  control  points 

3.  9509  unique  ground  points 

4.  Each  point  may  appear  on  up  to  10  photographs. 

5.  Object  space  is  expressed  in  space  rectangular 
coordinates  only. 

The  small  version  is  similar  with  the  dimensions  reduced 
to  150  photographs  and  150  ground  control  points.  The 
large  version  requires  476k  bytes  of  computer  memory  on  the 
IBM  370  computer.  The  core  capacity  is  reduced  to  282Kbytes 
in  the  smaller  version. 

Experience  has  shown  that  several  different  types  of  execu¬ 
tions  prove  helpful  in  data  analyzation.  The  idea  is 
similar  to  that  used  in  SIMBA.  One  must  try  to  separate 
the  affect  of  different  observations  in  order  to  locate 
blunders.  The  recommended  execution  sequence  is  as 
f ol lows : 


1.  the  i  n  L  o  r  s  o  c*  t  i  <>  n-on  I  v  option  with  no  '.round  -on- 

ttoi  .  Chock  tor  Lar-.i'  blunders. 

2  .  t;  s  c  the  tri.m"ulat  ion  option  with  no  ground  control. 
Chock  for  smaller  errors  in  point  mark  in;;  ind  measuring, 
especially  in  the  tie  points. 

S .  Use  tire  t r i an g u 1  a t i o n  option  with  "round  control. 
Weie.ii  the  image  measurements  tightly  and  allow  tiie  ground 
control  to  move.  Check  for  ground  control  errors  and  updat 
t  h  e  exposure  station  parameters. 

4.  Use  the  tri angulation  option  with  full  ground  control 
Weigh  the  solution  realistically;  check  all  residuals. 

5.  This  is  the  final  execution.  Use  the  triangulation 
option,  full  ground  control,  realistic  weights,  and  the 
latest  estimates  of  the  exposure  station  parameters.  Use 
the  error  propagation  option. 

There  are  several  things  to  watch  for  during  the  program 
executions.  Most  must  be  learned  through  experience.  One 
important  factor  to  remember  is  the  bandwidth  limitation. 

It  may  be  controlled  by  careful  ordering  of  the  frame 
position  and  attitude  cards.  One  can  stack  the  photos 
in  any  order.  The  trick  is  to  keep  the  frames  with  common 
points  as  close  together  as  possible.  The  pass  and  control 
point  selection  procedures  are  the  same  as  for  SIMBA  (see 
'figures  2  and  3).  It  should  be  noted  that  GIANT  will  not 
distort  block  corners  as  severely  as  SIMBA,  so  corner  posi¬ 
tioning  is  not  as  critical. 

The  output  listing  includes: 

1.  Camera  station  parameters  with  residuals. 

2.  Plate  coordinates  with  overall  standard  deviation  of 
x  and  y  per  f  rame . 

3.  Ground  control  coordinate  listing  with  residuals. 

4.  Camera  station  corrections  per  iteration. 

5.  Triangulated  residuals  per  frame  for  all  ground 
points. 

6.  Triangulated  camera  station  position  and  orientation 
residuals. 

7.  Triangulated  ground  point  coordinates. 

8.  Applied  ground  control  corrections. 

9.  Variance-covariance  matrix  and  standard  deviation  for 
the  camera  station  parameters  and  ground  control. 

Excerpts  from  a  typical  program  execution  may  be  found  in 
Append ix  F . 


As  wit!’.  SIMI’A,  problems  with  m  L  p  un  >•  h  e  d  or  m  i  jnunbo rud  data 
points  in  cl  errors  in  the  organization  o  i  the  data  deck  can 
cause  time  delays  and  needless  computer  expense.  A  data 
clean-up  program  was  written  by  the  author  in  1980  to  avoid 
such  problems.  This  program  (CTSRT)  is  similar  to  the  STMBA 
Sort  program  in  that  it  uses  the  exact  same  input  setup  as 
GIANT  and  does  a  point  by  point  analyzat  ion.  The  output 
consists  of  the  following: 

1 .  Data  deck  listing. 

2.  Frame  position  and  attitude  card  check. 

3.  A  frame  count  by  strip. 

4.  A  point  by  point  listing  to  include  all  locations 
encountered  and  bandwidth  computation. 

A  sample  of  the  program  output  may  be  found  in  Appendix  G . 

SUMMARY 

The  programs  presented  here  are  not  state-of-the-art  tech¬ 
nology.  They  do'not  have  the  most  rigorous  solutions 
nor  the  most  sophisticated  computer  manipulations.  What 
thev  do  have  is  a  good  combination  of  flexibility,  mathe¬ 
matical  rigor  and  computer  independence.  The  overall 
package  can  and  will  be  improved.  It  is  hoped  that  a 
new  adjustment  program  written  especially  for  mini-computers 
can  be  incorporated  into  the  system.  The  GIANT  program  is 
under  revision  to  make  it  more  user  oriented  and  to  add 
additional  options.  One  such  option  would  allow  output  in 
geographies.  The  package  is  presented  as  a  model  for  dis¬ 
cussion.  Its  greatest  testimonial  is  that  it  is  being 
extensively  used  and  has  proven  successful. 
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SAMPLE  OUTPUT 
FROM  CRAST 
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